Langevin Dynamics of Fluctuation Induced First Order Phase Transitions: self consistent Hartree 
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The Langevin dynamics of a system exhibiting a Fluctuation Induced First Order Phase Transition is solved 
within the self consistent Hartree Approximation. Competition between interactions at short and long length 
scales gives rise to spatial modulations in the order parameter, like stripes in 2d and lamellae in 3d. We show 
that when the time scale of observation is small compared with the time needed to the formation of modulated 
structures, the dynamics is dominated by a standard ferromagnetic contribution plus a correction term. However, 
once these structures are formed, the long time dynamics is no longer pure ferromagnetic. After a quench from 
a disordered state to low temperatures the system develops growing domains of stripes (lamellae). Due to the 
character of the transition, the paramagnetic phase is metastable at all finite temperatures, and the correlation 
length diverges only at T = 0. Consequently, the temperature is a relevant variable, for T > the system 
exhibits interrupted aging while for T = the system ages for all time scales. Furthermore, for all T, the 
exponent associated with the aging phenomena is independent of the dimension of the system. 



I. INTRODUCTION 



Type-II superconductors 1 1], doped Mott insulators |2], ultrathin magnetic films |3|,|4j,|5|], lipid monolayers |6], Raleigh-Benard 
convection |7], quantum Hall systems 1 8], are all systems that under appropriate conditions present stable phases characterized 
by the presence of modulated structures. The existence of these modulated structures is well understood on the basis of the 
Fluctuation Induced First Order Phase Transition Theory (FIFOT), first developed by Brazovskii |9[. This scenario predicts 
that systems in which the spectrum of fluctuations has a minimum in a shell in reciprocal space at a non zero wave vector, 
undergo a first order phase transition driven by fluctuations, in contrast to the second order transition predicted by mean field 
theory. Moreover, the strong degeneracy in the space of fluctuations induces the existence of many metastable structures at 
j— ( \ low temperatures, and since the experimentally observed structures are in general metastable, dynamical effects become very 
C — ■ ■ important. Unfortunately, the dynamical behavior of these systems is far from being understood. 

The existence and stability of metastable structures and the nature of the nucleation processes in the context of the Brazovskii 
scenario where first studied by Hohenberg and Swift |7']. They obtained the free energy barriers to nucleation and the shape 
and size of critical droplets in the weak coupling limit. Gross et al. lHoll compared the predictions of the self consistent Hartree 
approximation with direct simulations of the Langevin dynamics, confirming the validity of the approximation. 

A classic example where the Brazovskii scenario has got strong support, both theoretically and experimentally, is in diblock 
copolymers II llll 2lll 311 . These systems have interesting technological applications as self-assembling patterning media. Another 
*^ well known example of this kind of system is the three dimensional Coulomb frustrated ferromagnet. Wolynes et al. fl4ll have 
shown, using a replica dynamical mean field theory, that below a characteristic temperature, an exponential number of metastable 
states appears in the system preventing long-range order. Furthermore, through Monte Carlo simulations IU5II and also using a 
J-^ ' mode coupling analysis for the equilibrium Langevin dynamics of the Coulomb frustrated ferromagnet, Grousson et al. fl6j have 
J> . found an ergodicity breaking scenario in agreement with the predictions of Wolynes et al. These results resemble the behavior 
of many glass former systems, and two theoretical scenarios have been advocated in order to account for its phenomenology 
The relevance of the mode coupling predictions to the dynamics of the system have nevertheless been questioned by 
C$ • Geissler et al. O. 

The experimental and theoretical study of thin film magnetic materials have led to similar questions 1 20t l2ll l22ll . Thin films 
and quasi-two dimensional magnetic materials have many important technological applications, for example in data storage and 
magnetic sensors 12311 . In metallic uniaxial ferromagnetic films grown on a metallic substrate, like CoCu or FeCu, the system 
develops spontaneous stripe domains upon cooling below the Curie point. This is due to the competition between the exchange 
ferromagnetic short range interaction and the antiferromagnetic dipolar one, which is long range |4]. In the strong anisotropy or 
Ising limit, a self consistent Hartree approximation predicts the presence of a FIFOT for any value of the ratio between the ferro 
and antiferromagnetic coupling constants I21B . Nevertheless, the results from Monte Carlo simulations are far from conclusive, 
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suggesting the presence of first order transitions only for a restricted range of the ratio of the coupling constants J24ll25ll . The 
dynamics of these systems is also far from clear. Early work from Roland and Desai 12611 . who did simulations of the Langevin 
dynamics, concentrated on the early time regime, where modulation in the magnetization sets in. Later work with Monte Carlo 
simulations concentrated in some aspects of the out of equilibrium aging dynamics l27l Eill and the growing of stripe domains 
after a quench 1 29] . These works reveal a very rich phenomenology, with the appearance of complex phases reminiscent of 
liquid crystals, and strong metastability of the dynamics. All these facts point to the necessity of a more systematic study of the 
dynamical aspects of FIFOT. 

In this work we solve the Langevin dynamics of a generic model undergoing a first order phase transition driven by fluctu- 
ations. To characterize the long time dynamics of the system, we study the fluctuation spectrum close to the wave vector fco 
representative of the modulated phases. We solve the dynamical equations in the Hartree approximation and show that, already 
within this approximation, the dynamics of the system is very rich and departs from the usual ferromagnetic case. A key obser- 
vation is that, as we show in the next section, the spinodal of the high temperature disordered phase shifts to zero temperature in 
this approximation, and this has strong influence on the dynamics after a quench. In agreement with the equilibrium results, we 
show that the instability of the disordered phase appears only at T = 0, where the dynamics changes qualitatively. Nevertheless, 
the relaxation at finite temperature is far from trivial, showing the emergence of domains of stripes, which form a kind of mosaic 
state on top of the striped equilibrium phase. 

The rest of the paper is organized as follows. In section[I1]we present the model and show that, in the static self consistent 
approximation, it undergoes a FIFOT. In section[II]]we introduce the Langevin dynamics. In section HVl we present the general 
procedure to calculate the dynamical properties of the system in the Hartree approximation. Sections M and IVT1 are the core 
of the paper, where the results on correlations and responses are presented and discussed. Some conclusions are presented in 
section IVTI1 In two appendices we explain some technical details of the calculations. 



II. A MODEL WITH A FLUCTUATION INDUCED FIRST ORDER TRANSITION 

A classical model that undergoes a FIFOT may be defined by an attractive (ferromagnetic) short range interaction plus a 
competing, long-range repulsive (antiferromagnetic) interaction. In the simplest case of a scalar field, one can define an effective 
Landau-Ginzburg Hamiltonian of the form: 



■H[4>] 



d d x 



^(V0(f)) 2 + ^ 2 (f) + ^ 4 (f) 



1 

25 



d d xd d x' (f)(x)J{x,x')4>{x") 



(1) 



where r < and u > 0. J(x,x') = J(\x — x'\) represents a repulsive, isotropic, long range interaction and S measures 
the relative intensity between the attractive an d rep ulsive parts of the Hamiltonian. In the limit 6 — ► 00 one recovers the 
ferromagnetic O(N) model (for N = 1) OUHHHl- 

The (u/4) tp 4 term introduces a non-linearity that makes an exact solution of the model an impossible task. To deal with 
this non-linearity, one must consider the introduction of some kind of perturbative analysis. The simplest resummation scheme 
is the self-consistent Hartree approximation, or large N limit. It consists of replacing one factor of <fi 2 in the <fi 4 term in the 
Hamiltonian by its average (<j> 2 ), to be determined self-consistently. There are six ways of choosing the two factors of <fi to be 
paired in (</> 2 ), so the Hamiltonian in the Hartree approximation takes the gaussian form: 



H[4>] = t / d d x 



(V#?)) 2 +r^ 2 (f) + g<j> 2 (x) (0 2 (f)) 



where g — 3 it. Introducing the Fourier transform: 

4>(x) 
4>(k) 

the Hamiltonian takes the form 



d d k 



ik-x 



25 



4>{k) 



(2TT) d 

d d x e~ lts <j)(x) 



1 f d d k 1, * g f d d h 



d d k 2 



d d xd d x' 4>{x) J(\x- x'\) 4>{xi) 



(27r) d J (27r) rf J (2tt) 



d d k 



2 4>(ki)4>(k2)C(k 3 , -ki ~k 2 ~ k 3 ) 



1 

2 



d d k 



(2ir) d J (2tt) 
d d h f d d k 



~ d 0(h) 



Aik^S^^+g 



d d k 
(2-K) d 



C(k,-k-ki - k 2 ) 



4>{k 2 ) 



(27r) d J (2ir) d 



0(*l)A(*i,fc2)0(fc2) 



(2) 

(3) 
(4) 

(5) 
(6) 
(7) 
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with 

A(h,k 2 ) = A(h) 5 %x _ %2 + g J C (k, -k - fcx - k 2 ) (8) 

In the previous expression A(k) = r + k 2 + J(k)/8 and C(k, k') — (<f>(k)<f)(k')). Using that A(kx, k 2 ) = A(—ki, — k 2 ) and 
that C(k, k') = C c (k, k') + mp, where to^ = (<j>(k)) and C c (k, k') is the connected correlation function, we get finally the 
self-consistent Hartree equation for the connected correlator: 

A(k u k 2 ) = p^CzKhM = A(h) 5 %u _ %2 + g J {m %mii+ ^_ % + cSM+h-k)) (9) 

where (3—1 /ksT. In the paramagnetic phase, at high temperatures, all the order parameters = and the correlation matrix 
is diagonal, i.e. C c (k, k') — S c (k) _g, with S c (k) the static structure factor. From Eq.© we have 

P-'S-Hk) = A{k) +g Sc (k) = r + k 2 + ^±+g S c (k) (10) 

Introducing the "renormalized mass": 



(2n) d " K ' S 3 J (2tt) 



the structure factor becomes 



^ = T^Hf (12) 

where we have set kg = 1, and the renormalized mass has to be determined self-consistently from: 

/d d k 1 
-—7 T7u\ (13) 
(27T) d A + fc2 + ^l 

An instability in this equation may appear when A = A c = — (/c§ + j y ), where ko is the wave vector which minimizes A{k). 
Hence, the spinodal temperature T* is determined by the equation 



re-*-) j^m" k < 14) 



s 



where Kd is the surface of a d-dimensional sphere. The integrand in the right hand side is always positive and has a singularity 
at k = ko. Thus, the instability will be determined by the leading behavior of that integral, which can be estimated by expanding 
the denominator of the integrand around fco : 



f A fc"' 1 [ ko+e fc^ 1 r (k+koT^,, 

/ 77T7 « fc ~ / —f, r-^dk = / — dk = oo 

Jo X c + k 2 + ^ Jk -e c(fc-fco) 2 J-e k 2 

Therefore, the spinodal temperature always is depressed to zero, that is, (3*r* — > — oo. The fact that the isotropic phase is 
metastable at any finite temperature, a characteristic of the self-consistent nature of the fluctuations included in the model, will 
have important consequences on the dynamics after a quench at low temperatures. Nevertheless, it can be shown that below a 
melting temperature, the true equilibrium phase is a modulated one with characteristic wave vector ko and that a first order phase 
transition driven by fluctuations takes place @]. 



III. LANGEVIN DYNAMICS 



As usual, the Langevin dynamics for the scalar field <fi(x, t) is defined by: 



(15) 
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with, in addition, the following conditions for the thermal noise: 



(V(x,t)) = (16) 
(r)(x,t)r)(x',t')) = 2TS{x- x')5(t-t') 

In this work, we consider uncorrelated initial conditions: 

< rj(x, 0)77(2/, 0) >= A 8{x - x') (17) 

If 8 — ► 00, the last term in equation may be neglected and we keep only the short range part of the potential. At low 
temperatures this potential has two symmetric minima, 4>{x) — ±y / —^, where r < 0. The dynamics of this model is well 
understood|30, 31], and is determined, below T c , by the fixed point T = 0. At variance with the complete model, in the pure 
ferromagnetic case there is a continuous phase transition at a critical temperature T c which, in the Hartree approximation, is 
different from the mean field critical temperature, with a lower critical dimension d[ = 2. In this case, the dynamics after a 
sub-critical quench corresponds to an usual domain growth, with a growth law L(t) oc /j 1 / 2 , as in standard dynamical models 
with non-conserved order parameter. 

IV. SELF CONSISTENT HARTREE APPROXIMATION 
A. General Solution 

In our case of interest the dynamical equation ( TBI reads: 

d<j){x, t) = V 2^/^ t ) _ r ^/~ t) — u c/) 3 {x, t) - ^ f d d x' J(x,tf)<j>(tf,t) + r)($,t) (18) 
ot 2o J 

To extend consistently the known results for the equilibrium properties of the system to the study of the dynamics we keep the 
same resummation scheme. As before, in this approximation, the non-linear term <j> 3 is substituted by 3(0 2 (x, t))(f>(x, t) where 
the average is performed over the initial conditions and noise realizations. In such a way we obtain a linear equation in at the 
price of introducing a new parameter < <p 2 > to be determined self-consistently. To proceed, it is useful to go to Fourier space, 
in which we can write: 

^M = -[A(k)+I(t)]$(k,t) + fi(k,t) (19) 

where 

(»KM)) = (20) 

{fj(k,t)f)(k',t')) = 2T8(k + k')S(t-t') 



l(t)=r + g(0 2 (x,t)) (21) 



A(k) = k 2 + \j(k,P) (22) 




with initial conditions: 



(Mk)) = 
(Mk)Mk')} = (2ir) d A6(k + k') 



(23) 
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From equation ( 1191 it is easy to see that the general solution of the model may be written: 

(j>(k,t) = 4>(k,0) R(k,t,0) + ( R(k,t,t')fi(k,t')dt' (24) 

Jo 

where 

R(Ut') = ^e-^^ (25) 

andF(i) = e /o*'J(*'). 

Our main task is now to find a solution for Y(t), a function that encloses the unknown parameter introduced in the approxi- 
mation. Following standard procedures 1 30,^1 it is easy to show that: 

= 2rK(t) + 2gAf(t) + 4gT J dt'f(t - t')K(t') (26) 

where K(t) = Y 2 (t) and 

/(*) - / ^ 2ACk)t - < 27 > 

Equation J26i may be solve by Laplace transformation methods. If K(p) and f(p) are the Laplace transforms of K (t) and f(t) 
respectively, then equation (I26> reduces to: 

R = 2 fl A/-( P ) + g(0) 
p-2r-4gTf(p) 

Technically, the problem has been reduced to the calculation of /(p), to substitute it in equation (1281 and to calculate the 
corresponding Bromwich integral for K(t). Once K(t) is known, all the dynamical quantities of the system may be easily 
calculated from integral relations. 



V. CALCULATION OF K(t) 

In this, rather technical section, we go through a series of approximations and assumptions, which allow us to compute the 
function K(t) in the long time limit. We start the section presenting the approximation used to manage A(k). It keeps the 
necessary ingredients to model a fluctuation induced first order phase transition. We then proceed to the calculation of f(t) and 
finally K (t), in the long time regime. We show that the cases T = and T > give rise to different physics, in agreement with 
the static calculations. 



A. Approximation for A(k) 



Unfortunately, the analytical calculation of f(p) for general A(k) is a hopeless task. We will simplify it, considering only 
cases in which A(k) depends on the modulus of k, A(k) = A(k) (isotropic interactions). Since we are interested in the long 
time dynamics of the model, and we know that the equilibrium phases are characterized by the existence of a non-trivial wave 
vector &o ^ 0, at which the spectrum of fluctuations has a maximum, it is then natural to develop A(k) close to kg: 



A(k) =A Q + ^(k- k ) 2 + Q[(k - k y 



(29) 



where 



A = A(ko) 
d 2 A 
dk 2 

etc 



A 2 = 



k—kn 



(30) 
(31) 
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with A 2 > 0. 

Note that, if t is large enough, this approximation is valid not only for models with long range interactions, but in general is a 
good starting point to study other systems whose spectrum of fluctuations has an isotropic minimum at a non-zero wave vector. 
Therefore, the reader must keep in mind that the results of the next sections are valid in a context more general than the one 
represented by the Hamiltonian Q. 

From a technical point of view, one may note that Aq is irrelevant to the dynamical behavior of the system. In fact, from 
equation ( 119b . one can easily see that it is equivalent to a rescaling of r, and therefore to a shift in the critical temperature of the 
system. Therefore, from now on it will be neglected in our calculations. 



B. Results for /(p) 
With the previous assumptions for A(k), we may write f(t) as: 



/(*) 



27r d / 2 f d d k 



T(d/2) J (2ir) d 



-A 2 t(k-k ) 



whose Laplace transform becomes 



f(p) = 



d d k 



(2TT) d p + A 2 {k - k ) 2 



(32) 



(33) 



Next we analyze the behavior of this integral when p ~ 0. Adding a cutoff factor in the integrals in order to regularize the 
behavior for large wave vectors: 



m = 



d d k 



-(k-k ) 2 /A 2 



2n d / 2 



°° kd-l e -(k-k ) 2 /A 2 

dk 



(2ir) d p + A 2 (k-k ) 2 (2iT) d T(d/2) J p + A 2 (k~k a ) 2 
Then, after simple algebra equation (I34> becomes: 



(34) 



f(p) 



2n d / 2 



dke- pk ^ A2A ' 2 



ko 



_P_ 

A 2 



1/2 



(27r)^r(d/2) { v A 2 yi 2 J_^y" M 1 + k 2 

Expanding the binomial inside the integral we obtain the following expression for general dimensions: 



(35) 



I(P) = 



2n d / 2 



1 { A 2 \ 1 I 2 — 



(27r) d r(d/2) A 2 \p 



3=0 



d-1 

3 



3/2 

t 1 k 



d-l-j 



dke -pk 2 /A 2 A 2 



(#) 1/2 !M 



1 + fc 2 



(36) 



with d E Z. 

At this point two limit cases are possible and will be treated separately below. If (^■) 1 ^ 2 |fco | — ► the time scale of observation 
is such that the stripes are not completely formed, - < -hi-j-- In this limit we recover the dynamic properties of the pure 

P JT-2 

ferromagnet for fco ~ 0. On the other hand, if ( ) 1 fco I - * o° tne stripes are akeady formed and the interaction among them 
will be responsible of the dynamical properties of the system. Once this, more interesting limit, is taken, it is impossible to 
recover the pure ferromagnetic behavior. 



1. Stripes in formation 



Defining: 



J/2 

Fj(ko,p) =[±) h 



dke -pk 2 /A 2 A 2 

lTk 2 



k J 



(37) 
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and writing it as a Taylor series expansion, for fco ~ 0: Fj(ko,p) — Fj(0,p) + FU0,p)ko + O(k^), we get: 



f™ dkg^e-^ itj = d-l 







otherwise 



and 



JT dk&e-*b ifj=d-2 



otherwise 



Therefore, up to first order in k , F(ko,p) becomes: 



Fj(k ,p) = (^-) 



, oo 7 d—1 

- (JL\( d -v/ 2 1 - k — - 



o 1 + fc 2 



e A i a 



\4 



1 + Ar 



p fe 

e ^2 



and f(p) may be written as: 



(38) 



(39) 



(40) 



Hp) = V h 



fl + 6(^)i/2 fco 
V 



with a 



and 6 



c(¥) 



(41) 



(±n)HA 2 (i^Y^Al' 2 r(|) ' 

From equation J41i it comes out that the presence of modulated phases appear in the dynamics as a correction in f(p) to the 
usual ferromagnetic case [30]. One must remember, however, that this is true provided the second term within the brackets is 
small, i.e. during the formation of the modulated structures. 



2. Stripes formed 

On the other hand, for the limit (^■) 1 ^ 2 |fco| - * oo, stripes of sizes l/|feo| are already formed, and the dynamical properties 
of the system are defined by their interactions. The Taylor expansion, for fco finite and p — > 0, is written as: Fj(ko,p) — 
F J {k ,0) + F'(k ,0)p + O(p 2 ).Then, 



Fj(k ,0) 



h d-l roo dk _ ud-1 : f „■ _ n 

ft o J-oo I+F — nK o 11 3 — u 
otherwise 



and developing as before the first order derivative with respect to p, and taking the limit p — * 0, f(p) becomes: 



(42) 



f(p) 



1 (^ 2 -) 1 /2| 



(2ir) d r{d/2) A 2 y p 



TTk 4 ' 1 + A 2 k d ~ 3 



dk- 



1 



p k z 




(43) 



This last expression shows that for small p, f(p) — a + bp~ x / 2 with a < 0, independently of the system dimension. This limit 
was also explicitly calculated for d = 1, 2 and 3, confirming the series analysis. The calculations are shown in AppendixIXl 

Summarizing this subsection, f(p) was calculated in two limiting cases. In the first case, the system is still evolving and the 
stripes are not formed. The dynamical properties resemble the ones of the pure ferromagnet plus a correction term. Once the 
stripes are formed, the dynamic changes qualitatively, and one gets that f(p) = a + bp^ 1 / 2 , independently of the dimensionality. 
From now on, we will use this expression in future calculations, and only when necessary we will give explicit values for a and 
b. 
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C. The function Kit) 

By definition: 



1 />C7+200 

K{t) = — dpe**K(p). 

^ 7r * J a—ioa 

with K(p) defined by 1281 . where the regularization factors in f(p) can be disregarded in the long time limit. Then: 



(44) 



k( ) = 2gAa + 2gAbp~ 1 / 2 + l 

[P> p-2A ~2r-4gTa-4gTbp- 1 / 2 ( ' 

Because f(p) — a + bp~ x l 2 , K(p) has a branch point at p = 0, and the denominator varies in the domain (— oo, oo). 
Simplifying the notation we can write 

K( \ - A + Bp- 1 ' 2 Ap 1 ' 2 + B 

[P) p-C-Dp- 1 / 2 p z l' 2 - Cp 1 / 2 - D ' { ' 

with A = l + 2gAa, B = 2gAb, C = 2r + 4gTa and D = AgTb. 

The denominator of equation J46l > has three poles. Through a careful analysis it is possible to show that one pole is real and 
positive for all temperatures. The other two are complex conjugate with negative real part. We have to solve: 

1 r +!0 ° j_ a p 1/2 + b 



2«7 c _ loo * p3/ 2 ~Cp^ 2 ~D 

i r +io ° ... oPt a p 1/2 + B 



^j c _ loo dpeP ( P v 2 - x)( P i/ 2 - z)( P i/ 2 - z *y (47) 

where x 2 £ 1Z, z 2 £ C and (z 2 )* is the complex conjugate of z 2 . After a lengthy computation (see AppendixlEl. we get: 

= 2x(Ax + B) ^ H 2z(Az + B) 2 z* z* + g) ^ yt 

(x — z)(x — z*) (z — x)(z — z*) (z* — x)(z* — z) 

IT Br3/=» + (gC->U)r^ 

7ri ^2 + (r 3/2 +Cr l/2)2 ^ 

Now, one must distinguish carefully the cases T > and T = 0. Note, for example, that the real parts of the complex 
poles are negative even for T — > 0, while the real pole is always positive going to zero at T = 0, where the physics changes 
qualitatively. 



1. T > 

In this case one may neglect the contributions from the complex conjugate poles, since their real parts have decaying expo- 
nential functions. On the other hand, the last integral in equation (I48> may be easily estimated noting that: (r 3 / 2 + C r 1 / 2 ) 2 = 
r 3 + 2 C r 2 + C 2 r, and that for long times t — > oo the dominant contributions will come from r -C 1. We end with: 



AC-^ r dr e -H r i/2 _ BC-AD ^ 
D 2 J D 2 2t 3 / 2 



(49) 



Therefore, 



1 ' ~ (a: - z) (z - z* ) 2 V^D 2 i 3 / 2 ' P ' 



where the last term goes to zero for t — > oo. 
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2. T = 



In this case one must note that x(T) — ► T, therefore the contribution from the real pole dissapears. At the same time, the 
complex poles converge to a single real pole that gives rise to a decaying exponential function. 

Moreover, the expansion used to calculate the last integral in ( I48> is not longer valid. Being D = one finds that, for large t: 

roc - rt B 

B / dr 3 ~ — — - (51) 

and 

„ rt b r(i/2) 

with C < 0. The first term comes from the limit as T — > of the two complex poles and is obviously subdominant in this 
analysis. 

Already with these results at hand one must note two important differences with the usual ferromagnetic coarsening. The 
first one is that here the temperature is a relevant variable. While for T > the relaxation will be dominated by exponential 
(paramagnetic) contributions, for T = the relaxation will be power-like. The second one is that, excluding irrelevant prefactors, 
the long time dynamics, for all temperatures, is independent of the system dimensionality. 



VI. RESPONSE AND CORRELATION FUNCTIONS 



In this section we present the main physical results of the paper, regarding the behavior of correlation and response functions. 
As seen in Section dIV> . the two-times response function in Fourier space is given by: 

fl( M ,i') = ^ e -^)(*-*') (53) 



with Y(t) = y/ K{t). Defining the two-times structure factor: 

{4>{k, t) <j>{k',t')) = (27r) d S(k + k') C(k, t, t') (54) 
and using i20\ . (12 1 i and 11241 . one can show that: 

C(k, t, t') = A R(k, t, 0) R(k, t', 0) + 2T f R(k,t,s)R(k,t',s)ds (55) 

Jo 

As the physics at finite temperature is different from that at T = 0, we will analyze both cases separately. Also, note that, for 
T ^ 0, the leading contribution to K(t) is exponential, and consequently, to leading order, correlations and responses will be 
stationary, consistent with a paramagnetic phase. Nevertheless, the algebraic subdominant contribution precludes the presence of 
a transient non-stationary dynamics, with time scales that can be large for low enough temperatures. Consequently, we will keep 
the subleading contribution also and analyze its effect on the dynamics, showing that it leads to interrupted aging in correlations 
and responses. 



A. T > 

From equation J50i we find: 



{ ) y (x-z)(x-z*) 2V^D 2 i 3 / 2 ' 

and defining C 2 = ^-z){oc^) Ci = — ^^_pf to simplify the notation, one gets for large t: 



Y {t)=cle^{l + ^^) (57) 
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Then, the response function for T > is: 



R(k,t,t') = e -hl* 2 +M(k-k a f\(t-t') 



1 + 



C 2 e- xH ' 



(58) 



2Ci t' 3 / 2 

The two-time correlation function is given by equation (I55i . Then, using fl58i and defining B(k) — x 2 + A 2 (k — fcp) 2 we get 



C(k,t,t') = Ae-5 fl W(*+*') + 2Te- = - B W(*+*') / ds e~i B ^ s 1 + 

Jo \ 2C i s 7 

Performing the integration in ( I59> and setting t = t — t' one gets: 



(59) 



C(k,r,t') = e ~iB(k)T 



2T 



A 



2T 
BQcj 



-B(k)t' 



(60) 



The second term within parenthesis goes rapidly to zero, while the \/V in the third term reflects the presence of interrupted 
aging in the system. On the other hand, keeping the stationary part and for t = t', we obtain the static structure factor: 



2T 2T 
C{k) = lim C(k, t) = — - = 



(61) 



which, as expected, shows a characteristic peak at k = k$. 
The correlation function in real space is given by: 



f°° d d k 

C W= Trh c ( k )e ik ' x ( 62 ) 



In d dimensions: 



d-2 



1 f°° I 1 \ k d ~ l dk 

c w = jw&J [vr) J ^ {kr) x 2 + A 2 (k- ko r (63) 



where J„(x) is a Bessel function of the first kind. In the limit kr — > oo, 



d-2 

,<i-l . 



C{r)K j2^J {Vr) V^ C0S r"~4~J ^+A 2 (fc-fc ) 2 (64) 
This integral can be solved using the theorem of residues in the complex plane. The final result is: 

e -r/t 

C(r) ex cos {kgr — tp) d _ 1 . (65) 

7* ~ 

where ip = ( d ~ 1 ^ ,r — tan -1 ( 9 ^~Vf ) • We see the presence of a correlation length and a modulation length. The correlation 



4 V 2feo v^i2 

length is given by 




As mentioned before, at low temperatures, x(T) oc T, and consequently the correlation length diverges at T — as £(T) oc 
1/T. Then, one can conclude that after a quench from the disordered phase to a very low temperature the system breaks into 
regions inside which there is modulated order, or stripes. No long range order is observed. A transition is approached at T = 0, 
where the correlation length diverges and stripe order sets in. In this respect, it is interesting to see also what happens with the 
correlation function at T = 0. 
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B. T = 



For T = we may neglect the decaying exponential and simple algebra gives: 




(67) 



Substituting d67i in (I55> it is easy to prove that 



C(k,t,t') = AR(k,t,0)R(k,t',0) = A(U')ie-^^- ko ^^ t+t '^ 



(68) 



and making t = if we obtain 




-A 2 (fc-fc ) 2 (t) 



(69) 




For d = 1 the spatial correlation function has the form: 



C(r,t) 



cos (fegr) e 



4t g r 



.2 



f(k Vt), 



(70) 



where erf(x) is the error function. This result is consistent, for large t, with the appearance of modulated structures and 
long-range order. Nevertheless note that, from Ea. (l65> . a critical decay of correlations is observed for d > 1. 

The stability of the paramagnetic phase until T — is obtained also in a purely static calculation of the phase diagram in the 
self consistent (Hartree) approximation. Nevertheless, the same calculation shows that the disordered phase is only metastable, 
the stripe phase has a lower free energy below a finite critical temperature and is thus the true thermodynamic equilibrium of 
the system. Our calculations show that the Langevin Dynamics within the Hartree approximation reproduces this scenario. A 
quench from a disordered phase to T > gives rise to a paramagnetic-like dynamics reflecting the metastability of this phase. 
Only at T = 0, the spinodal is reached, and the system shows a coarsening non-equilibrium dynamics, with diverging time 
scales in the thermodynamic limit. 



In this work we have calculated within the Hartree approximation the exact long time dynamics of a model system exhibit- 
ing a Fluctuation Induced First Order Phase Transition. We motivate our work starting with a Hamiltonian with short range 
ferromagnetic interaction and long range antiferromagnetic interactions, but our results are valid in general for the long time 
dynamics of any system exhibiting Fluctuation Induced First Order Phase Transitions, provided that the spectrum of the fluctu- 
ations is isotropic. We present explicit expressions for one and two time correlation and response functions and show that the 
dynamics converges to known static results in the Hartree approximation. Our results show that the dynamics of the system may 
be decomposed in two stages, first the modulated phases form, and during this stage the dynamics follows the usual coarsening 
scenario for a ferromagnetic system: it is dominated by the zero temperature fixed point and depends on the dimensionality of 
the system. Once these modulated structures are formed, the dynamics changes qualitatively. It becomes independent of the 
system dimension and the temperature becomes a relevant variable. For T > the system exhibits interrupted aging and a 
standard paramagnetic relaxation for large times, dominated by the presence of metastable states. At low temperatures domanins 
of stripes are formed. At T = the correlation length diverges and stripe order sets in. The system ages for all the time 
scales following a coarsening dynamics that searches the equilibrium state. Moreover, the exponent associated with this aging 
dynamics is independent of the system dimensionality. 

In this work we have explored only the presence of positional order, through the calculation of the correlations of the field 
4>(x,t). It would be interesting to compute also orientational observables, which are known to be relevant for these kind of 
systems, and give rise to nematic-like order |25]. Other interesting questions that can be addressed starting from the present 
calculations are the possible nucleation of stripe phases in the paramagnetic state, which eventually should lead to the first order 
transition predicted within the static Hartree approximation. Also the possible presence of freezing in the low temperature dy- 
namics could be addressed withing a refined approximation, like mode-coupling or the self-consistent screening approximation. 



VII. CONCLUSIONS 
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APPENDIX A 



Here we show the explicit calculation of f(p) in the limit p — > for d = 1, d = 2 and d = 3. We begin from Eq. (I36> in each 
case, and develop the sums and integrals. 



d=l 



f(„) = i r dk e -v^iA^ (A1) 



Now, 



^-pk 2 /A 2 A 2 = 



= 7T- / _T_ e -Pfe 2 M2A^ (A2) 



The last equality is due to the integrand be an even function. For p — > the last integral gives: 

(A3) 



/ dfc pfc 2 M2 A 2 = f A ] £ + Of„3/ 2) 



Then, 



1 f / v \ 1/2 e-l fe °l 2 /A 2 , 1 

i -|& | 2 /A 2 

- ' |, | + 0(p) (A4) 



4/2" 7rA 2 |fco 



d=2 



/(P) 



r/A: 



2 7 r(M 2 ) 1 / 2 7-(^) 1/2 | fco | 1 



-pk 2 /A 2 A 2 



k 2 



p \ 1 ^ 2 k 



A 2 



Ao 



Here we have to solve two integrals. The first one was already solved for the case d = 1: 



dk 



(4?-) 1/2 im 1 + fc2 



e -pfc 2 /A 2 A 2 = ^ 



i) T +ofr/) 



The second integral is: 



(A5) 



(A6) 



dk k 



//l2\l/2|, 1 



-pfc 2 /A 2 A 2 _ 



dfc A 



-(^) 1/2 |fc | 1 + 7_oo ./-oo J^ 2 lkol l + k 2 

because the integrand is an odd function. The last integral can be approximated as: 



-pfc 2 /A 2 A 2 



dkk 



-pfc 2 /A 2 A 2 _ 



dA 



-pk 2 /A 2 A 2 



(4^) 1/2 im 



(A7) 



(A8) 
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Now 



(#) 1/2 !M k 



^ e -pk 2 /A 2 A 2 



(A9) 



where 7 is the Euler constant. The other integral 



Finally, 



(#) 1/2 |M fc3 



* - P fc 2 A4 2 A 2 



1 r pe -iM 2 /A 2 p r(o,l$-) 



A 2 |fc | s 



Ifcnl 2 
A 2 

^ 2 A2 



P I £ 
2A 2 I 



l*b| 



(A10) 



Then, 



/oo 



-pfc 2 /A 2 A 2 



7 1, (W 2 
■log 



2 2 fa V A 2 / 2A 2 N 2 



p e 



-|feo| 2 /A 2 



+ 0(A- 



(All) 



/» 



fc 



27 rA 2 /2 



Tip 



-1/2 



X e -|fc | 2 /A 2 x 
+ 



4 /2 1*0 



k 



2A 



1/2 J 



-1/2 _ 



2ttA 2 



AttA- 



-OQogA- 



(A12) 



d=3 



In this case 



2tt 3 / 2 A; 2 



(2^)3r(3/2)(pA 2 ) 1 / 2 7-(^) 1 / 2 | fc0 | 1 + fc 2 



-pk 2 /A 2 A 2 



27T 3 / 2 fc 2 



(2^)3r(3/2)(pA 2 )V2 



_P_ 

A 2 



1/2 e -|fe„| 2 /A 2 j 



1*0 I 



+ 



1 + 



p 



ko \A 2 



P_ 

A 2 

1/2 



1/2 



k 



p_\ fc 2 
A 2 ) k% 



7 1 , 

log 

2 2 5 



\kol 
A 2 



fr 2 



A 2 



dk fc 2 

1/2 i fc „i TTfc 2 "' 



-pk 2 /A 2 A 2 



(#) ' iM 



(A13) 



Again, the contribution of the dominant term is of order p 1 I 2 , proving that, in agreement with the series expansion in the 
text, the behavior of f(p) is independent of dimensionality. 



APPENDIX B 

For f(p) == a + bp~ x l 2 we have 



1 [ c+i °° pt Ap 1 / 2 + B 



K(t) = — dpe 



1 f C+l °° ApV2 + B 



= 2^7™ dpe " (pVa-aOfri/a-*)^-^)' (B1) 
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FIG. 1: Path of integration for the inverse Laplace transform 



where x 2 E 1Z, z 2 6 C and (z*) 2 is the complex conjugate of z 2 . The three poles x 2 , z 2 and (z*) 2 are simple and the residues 
are: 

hm (p - x 2 ) e pt K(p) = ±- -f — 

p^x 2 [X — Z)\X — Z ) 

lim (p - z 2 ) K(p) = ( f Z + ^ eZ2t * Z (B2) 
p-^z 2 (z — x)(z — z*) 



lim (p- (z*) 2 ) e pt K(p) 



(Az* + B) e( z *) 2 * 2z* 



p^{z*) 2K (z* - x){z* - z) 

In order to perform the integral, consider the contour in Fig.Q. We choose the branch cut to be the negative real axis (— oo, 0]. 
From JB 1> we conclude that the paths BC, DEF and GA do not contribute. Then: 



t'B t>D t>G 

/ + / + / = 2m V Res 

J A JC JF 



(B3) 



Define p — re 1 



° , nt Ap^ + B f\. _ pe «„ t i^e^ + B 



r p3/ 2 — Cp 1 / 2 — D 7^ V ' r 3/2 e ai 1 r/2_ C ' r l/2 e iV2_i5 



° dre - * + ^ (B4) 



p3/2 _ (7^1/2 _ ]j Jp y ' r 3 / 2 e - 3i7T / 2 - C r 1 / 2 e - i7T / 2 - D 

B — iAr 1 ! 2 

dTe ^ ir 3 / 2 +iCr^ 2 -D (B5) 



o 
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Then, 

' ;; ' x ( -B + iAr 1 / 2 B + iAr 1 ' 2 



JC JF JO I 



-I 

-i: 



oo 

rt 



dr e 



i(r 3 / 2 + C r 1 / 2 ) — D i(r 3 / 2 + C r 1 / 2 ) + D 

B-iAr 1 ' 2 B + iAr 1 ' 2 

D - i(r 3 / 2 + C r 1 / 2 ) ~ D + i(r 3 / 2 + C r 1 / 2 ) 
[D + i(r 3 ' 2 + Cr 1 ' 2 )]{B - lAr 1 ' 2 ) - [D - i(r 3 / 2 + Cr 1 l 2 )]{B + lAr 1 / 2 ) 



D 2 + (r 3 / 2 + Cr 1 / 2 ) 2 
The numerator has the form: 



(B6) 



xy*-x*y = (Mx + i%x){My - iSy) - (Mx - i$sx)(?lty + i$sy) 

= 5i:z:3fy - mxSy + mySx + SxSy - iVixZiy + m y %x - $ix$iy - ^sx^sy 

= -2mx^y + 2my^sx (B7) 



Then, 



J + I 



F JO 



7 



D 2 + (r 3 / 2 + C r 1 / 2 ) 2 
Br 3 / 2 + (BC - ADy' 2 



2i dre~ rt -J— , (B8) 



D 2 + (r 3 / 2 + Cr 1 /2) 



and 



ACKNOWLEDGMENTS 

We gratefully aknowledge partial financial support from the Abdus Salam ICTP through grant Net-61, Latinamerican Network 
on Slow Dynamics in Complex Systems. D.A.S. thanks Masayuki Hase for helpful discussions and CNPq, Brazil, for financial 
support. 



[1] J. M. Tranquada, B. J. Sternlieb, J. D. Axe, Y. Nakamura, and S. Uchida, Nature 375, 561 (1995). 

[2] S. A. Kivelson, E. Fradkin, and V. J. Emery, Nature 393, 550 (1998). 

[3] O. Portmann, A. Vaterlaus, and D. Pescia, Nature 422, 701 (2003). 

[4] K. De'Bell, A. B. Maclsaac, and J. P. Whitehead, Rev. Mod. Phys. 72, 225 (2000). 

[5] E. A. Jagla, Phys. Rev. E 70, 046204 (2004). 

[6] M. Seul and D. Andelman, Science 267, 476 (1995). 

[7] P. C. Hohenberg and J. B. Swift, Phys. Rev. E 52, 1828 (1995). 

[8] D. G. Barci, E. Fradkin, S. A. Kivelson, and V. Oganesyan, Phys. Rev. B 65, 245319 (2002). 

[9] S. A. Brazovskii, Sov. Phys. JETP 41, 85 (1975). 
[10] N. A. Gross, M. Ignatiev, and B. Chakraborty, Phys. Rev. E 62, 61 16 (2000). 
[11] G. H. Fredrickson and K. Binder, J. Chem. Phys. 91, 7265 (1989). 
[12] F. S. Bates, J. H. Rosedale, and G. H. Fredrickson, J. Chem. Phys. 92, 6255 (1990). 

[13] D. E. Angelescu, C. K. Harrison, M. L. Trawick, R. A. Register, and P. M. Chaikin, Phys. Rev. Lett. 95, 025702 (2005). 

[14] H. Westfahl Jr, J. Schmalian, and P. G. Wolynes, Phys. Rev. B 64, 174203 (2001). 

[15] M. Grousson, G. Tarjus, and P. Viot, Phys. Rev. E 64, 036109 (2001). 

[16] M. Grousson, V. Krakoviack, and G. Tarjus, Phys. Rev. E 66, 26126 (2002). 

[17] J. D. Stevenson and P. G. Wolynes, preprint cond-mat/0505070 (2005). 

[18] G. Tarjus, S. A. Kivelson, Z. Nussinov, and P. Viot, preprint cond-mat/0509127 (2005). 

[19] P. L. Geissler and D. R. Reichman, Phys. Rev. E 69, 021501 (2004). 

[20] T. Garel and S. Doniach, Phys. Rev. B 26, 325 (1982). 

[21] S. Cannas, D. Stariolo, and F. Tamarit, Phys. Rev. B 69, 092409 (2004). 

[22] O. Portmann, A. Vaterlaus, and D. Pescia, Phys. Rev. Lett. 96, 047212 (2006). 



16 



[23] A. Hubert and R. Schafer, Magnetic Domains (Springer- Verlag, Berlin, 1998). 
[24] E. Rastelli, S. Regina, and A. Tassi, Phys. Rev. B 73, 144418 (2006). 

[25] S. A. Cannas, M. Michelon, D. A. Stariolo, and F. A. Tamarit, Phys. Rev. B 73, 184425 (2006). 
[26] C. Roland and R. Desai, Phys. Rev. B 42, 6658 (1990). 

[27] J. H. Toloza, F. A. Tamarit, and S. A. Cannas, Phys. Rev. B 58, R8885 (1998). 
[28] D. A. Stariolo and S. A. Cannas, Phys. Rev. B 60, 3013 (1999). 

[29] P. M. Gleiser, F. A. Tamarit, S. A. Cannas, and M. A. Montemurro, Phys. Rev. B 68, 134401 (2003). 

[30] F. Corberi, E. Lippiello, and M. Zannetti, Phys. Rev. E 65, 046136 (2002). 

[31] T. J. Newman and A. Bray, J. Phys. A: Math. Gen. 23, 4491 (1990). 

[32] C. Chamon, L. Cugliandolo, and H. Yoshino, preprint cond-mat/0506297 (2005). 



